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1 . 0 INTRODUCTION 

The purpose of this study is to document the sensitivity of 
multivariate statistical interpolation to the values of 
statistical parameters, misspecif ication of these parameters, and 
to missing observations. It is hoped that this information will 
enable practitioners to concentrate on the appropriate 
specification of crucial parameters while avoiding agony over 
those that have relatively small effect. The results presented 
here point out some interesting dependencies which I believe have 
not been previously published. Prior studies upon which this 
work builds are Franke (1985 and 1988) and especially Franke, et 
al . (1988). It was also influenced by Seaman (1983). These 
previous studies, however, involved only univariate objective 
analysis schemes. While many of the results can be expected to 
carry over in a similar way, the influence of wind observations 
on the expected error over various observation sets is useful and 
interesting. Further, I believe that the graphical presentation 
of the results used here makes it easy to discern the important 
parameters . 

A brief overview of multivariate statistical interpolation 
and the method used to calculate expected errors is given in 
Section 2. The details of the necessary calculations are not 
given explicitly there, but are available in Appendix 2, which 
gives a listing of the subroutine used to evaluate the covariance 
matrices. Section 3 gives a discussion of the various parameters 
and the values over which they were varied during the study. 
Section 4 contains an analysis of the results. Because of the 
plethora of data which was generated it is difficult to 
comprehend the important details in tabular format, and therefore 
the information is incorporated into a few graphs which enable 
one to easily ascertain the sensitivities of the scheme to the 
various parameters. Tables are included in Appendix 1 for 
completeness, but it is expected that few readers will find them 
necessary. Section 5 contains the concluding remarks. 
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2.0 MULTIVARIATE STATISTICAL INTERPOLATION 

Statistical Interpolation (SI) is in use at many of the 
world's Numerical Weather Prediction (NWP) facilities, including 
the U.S. Navy's Fleet Numerical Oceanography Center (FNOC). The 
scheme has its roots in the work of Kolmogorov and Weiner, and 
was first developed for meteorological applications by Gandin 
(1965). More recently it has been applied in multivariate form 
by Schlatter (1975), Schlatter, et al. (1976), Bergman (1979), 
and Lorenc (1981). 

In its perfect form (all parameters known), SI delivers 
estimates of a field with minimum mean squared error over a 
certain ensemble of realizations that satisfy normality and 
stationarity of the underlying stochastic process. For 
convenience, isotropy is usually assumed, and for meteorological 
problems a zero mean is assumed (although a nonzero mean can be 
accounted for in more than one way) . 

In the multivariate formulation, the dependent data consists 
of related variables, which in our case we assume to be the 
errors in the background pressure height and wind fields, H, U, 
and V, as related through the geostrophic relationship: 

U - V - kjHy . 

Here the values of kj^ and are dependent on the latitude. The 
SI equations are applied to the background error, obtained by 
forming the difference between the background (normally the NWP 
values, interpolated to the observation locations) and the 
observed values. The assumption of normality implies the minimum 
variance (or least mean squared error) predictor is a linear 
combination of the data values. Then construction of the weights 
leads to solution of a system of equations whose coefficient 
matrix is the matrix of spatial covariances and cross-covariances 
for the variables. Because of the assumed relationship between 
the variables, the wind error covariances and the cross- 
covariances are determined when the height error spatial 
covariance function is known. The exact relationship is given in 
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Franke, et al . (1988), as well as in other references above. 

The SI equations have been derived in numerous places (see 

e.g. , Schlatter (1975), Lorenc (1981), Thiebaux (1985), and 

Thiebaux and Pedder (1987)) and are repeated here only for 

completeness. Let O = (o.) represent N measured values, 

J J~-*- / • • ’ / w 

with corresponding independent observation error variance 
. . 2 

(diagonal) matrix S = {a^}, B = (b^) the corresponding vector of 

background values, and C = (c. .). „ the spatial 

covariance matrix for the background errors. Then, letting 

C,^={c_.), . , „ denote the vector of covariances between the 

0 0 3 3=1 , . . . , N 

error in the variable at the location at which it is to be 
analyzed and the background errors, the following equation holds 
for the analyzed value, a^: 

a^ = b^ + cJ(C+E)"^(0 - B) . 

If the statistical parameters are known precisely (whereupon 
SI becomes Optimum Interpolation, or 01), then the expected error 
variance for the estimated variable at P^ is given in the usual 
least squares form, 

2 

where is the variance of the background error. Considering 
the case where the statistical parameters are not known exactly, 
the analyzed value becomes 

a^ = bQ + cJ(C+S)"^(0 - B) , 

where the tilde overbar signifies assumed inexact values. Notice 
that the equation is exactly the same except that assumed values 
replace the exact values for covariances. The expected mean 
squared error is 

- 2cJ(C+S)"^Cq + cJ(C+E)"^(C+E)(C+E)"^Cq . 
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These equations are given in Seaman (1983), where the vector 
Wq = (C+E ) is interpreted in the usual statistical 

fashion, as the weight vector for the observed minus background 
values. This also simplifies the writing of the error expression 
and indicates a computationally efficient algorithm. 

The expected error will be computed at a number of 

locations, in our case on a 7x11 grid of points. The matrix 

(C+E) is symmetric and positive definite, so the computation of 

W is accomplished by performing a one time Cholesky 

~ . T 

decomposition of (C+S) into LL , followed by forward and 
backward substitutions to find W as the solution of 
LL Wq = Cq . This IS an important concept since the number of 
equations is the number of observations (up to 108 here), and 
varies over 77 locations in our case. In actual practice Pq 
would vary over several locations for the same observation set 
when a block analysis scheme is used. 



3.0 SETTING FOR THE STUDY 

An empirical study such as this one requires a number of 
compromises concerning the range of parameters permitted. In 
previous studies (Franke, 1985, and Franke, et al . , 1988) a set 
of three grids and corresponding observation locations based 
approximately on the radiosonde network over the United States 
and the Atlantic Ocean were used. The study here is also at a 
single level. One difference from my previous studies, which has 
a rather minor influence, is that distances are measured in 
meters (approximately; I did not use the exact geodesic 
distance) rather than degrees. The formula used for the distance 
(also used by Schlatter (1975)) between two points (0 = 
longitude, tp = latitude) is 

2 2 2 2 2 *^-i~*^i 

= r'"[(.pj-(p. )^ + (0j-e. )^cos^(-l^)] . 

While computed distances are (relatively) quite different in the 
upper latitudes, the overall influence on the expected mean 
squared errors is rather small over the regions we consider. 
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Figures 1-3 show the 2.5° grids and observation locations 
used in this study. The +'s indicate the grid points, while the 
small squares and circles represent the observation locations. 

The shaded circle represents the "missing observation" for tests 
of sensitivity of the analyzed values to one missing observation 
(both missing wind only and missing height and wind). The open 
circles represent additional missing observation locations for 
tests of sensitivity to many (one-half) missing observations. 

Each grid is a 7x11 grid of points, taken to be the interior grid 
points of a 9x13 grid containing the observation locations. The 
Middle United States (MUS) grid has 36 observations, the East 
Coast (EC) grid has 25 observation locations, while Middle 
Atlantic (MA) grid has only 3 observations. 

The assumed correlation function for the height errors is 
the specialized second order autoregressive function (SOAR) that 
seems to be quite stable with respect to variations in the 
parameters while embodying enough parameters to fit historical 
data (Franke, et al . , 1988). As noted by Thi6baux, et al. 

(1986), Balgovind, et al. (1983), the SOAR seems to have an 
affinity for meteorological data. The form we used is 

OT 

C(s) = (l-A)(l+as)e + A, 

where s denotes the distance, and a and A are parameters which in 
practice are determined by a fitting process. In the previous 
section c. .=C(d. •)af, where d. . is the distance between the i 
and j observation points. The study conducted here made use of 
five different sets of parameter values, as shown in Table 1, and 
also depicted in Figure 4. The "nominal" correlation function is 
considered to be #4, which approximates the Bessel function curve 
of Lonnberg (1982) closely, and also corresponds closely to the 
decay rate used by Lorenc and Mammon (1988). Varying from one 
correlation function to the next makes it possible to determine 
whether the additive constant or the decay rate constant is the 
more critical. 
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Cor# 


A a 1 


1 #m to #m+l 


1 


0.0 5.000X10“^ 




2 


0.0 3.082xl0~® 


a decreases 


3 


0.15 3.082X10“^ 


A increases 


4 


0.2722 3.082x10“^ 


A increases 


5 


0.2722 2.188x10“® 


a decreases 



Table 1: Correlation Functions 

The nominal value for the standard deviation of the 
background height error is 30 m. The correlation function then 
determines the variance of the wind field errors, and the 
standard deviations for wind errors are given in Table 2 . 

Because the values of and depend on the latitude, the rms 
errors over a given grid depend on the grid. All expected wind 
errors were calculated relative to their value at a given point, 
and the rms values of these are given in the tables and the 
figures, different than computing the expected rms values over 
the grid and then comparing this with the rms background error 
over the grid. 

The nominal values assumed for the observation errors were 
10 m for heights and 1.0 m/sec for wind components. For 01 these 
values were varied over the values 0, 5, 10, 20, and « m for 
heights, and 0, 0.5, 1.0, 2.0, and «> m/sec for the winds. The 
expected rms error was not computed for all combinations of these 
values; the « values imply no measurement of that variable and 
computations were only performed with the nominal value of the 
other measurement error, for example. These calculations were 
performed with 01 to show the effect of no observations of a 
particular variable. With no wind observations, the process 
collapses to the usual univariate 01 scheme. 

In order to assess the effect of missing observations, four 
"missing observation" computations were performed. For one 
observation point in each grid it was assumed that the wind 
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observations were not available; then it was assumed that the 
height observation at the same location was not available. For 
about one-half of the observation points it was assumed that the 
wind observations were not available; then it was assumed that 
the corresponding height observations were not available. 

The computations of expected rms errors for the statistical 
schemes did not cover such a wide range of assumed observation 
error values. Variations of the assumed observation errors and 
spatial correlation function number were only carried out to the 
adjacent value (assumed observation error changing by a factor 
of .5 or 2, and correlation function number only to the adjacent 
number) . 



Cor# 1 


1 grid 1&2 


grid 3 


1 


16.05 


14 . 52 


2 


9.09 


8.95 


3 


9.12 


8.25 


4 


8.44 


7.64 


5 


5.99 


5.42 



Table 2; Background wind errors, m/sec 
4.0 ANALYSIS OF THE RESULTS 

The easy detection of the sensitivity to parameters on which 
a process depends is sometimes clouded by the mass of information 
available. In the following, I believe the plotting scheme 
adopted enables readers to easily detect critical parameters. 

In this study there are four important parameters that are 
varied for each of the three grids: The height observation 

accuracy, the wind observation accuracy, the additive constant A 
in the correlation function, and the decay rate a (the reciprocal 
can be thought of as some measure of "correlation distance") in 
the correlation function. In this study the correlation function 
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parameters have not been varied except one at a time, in the 
manner noted in the previous section . 

We first consider the sensitivity of 01 to the accuracy of 
the observations. The results for the error in the analyzed 
height field for 13 different combinations of observation 
accuracies and 5 correlation functions are shown in Figures 5-7 
for the MUS, EC, and MA grids, respectively. Note that only the 
integer abscissa values have meaning, and the points are 
connected only to enable one ro more easily see the effects of 
changing parameter values. As one moves to the right on the 
scale, the trend is toward less observation accuracy, with wind 
observation accuracies first decreasing, while the height 
observation accuracy more slowly decreases. From the plots, it 
is seen that significant increases in 01 height analysis errors 
occur at abscissae 2, 5, 9, and 13, while between there is 
generally some increase, but much smaller in magnitude. Table 3 
lists the values of the observation accuracies, and it is seen 
that for 2, 5, 10, and 13 there are jumps in the observation 
accuracy for heights. The jump in expected error at abscissa 9 
is due to the complete loss wind observations, and if this 
abscissa is omitted, then the jump would be due to the height 
accuracy change at abscissa 10. From this graph it is apparent 
that the accuracy of height observations are of premier 
importance in the height analysis, while the accuracy of the wind 
observations are less important in the nominal range considered 
here, except that not having wind observations at all also 
results in a significant increase in 01 error. 

Significant increases in error of the analyzed height field 
are seen to occur with decreasing correlation distance 
(increasing value of a) and decreasing constant A, for all three 
grid and observation point sets. Of course, the effects of 
observation accuracy are considerably smaller for the sparse MA 
grid. For all three grids it is apparent that 01 errors are more 
sensitive to changes in the correlation distance (parameter a) 
than in the additive constant A. 
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abscissa 


h error 


w error 


1 


0 


1.0 


2 


5 


0.5 


3 


5 


1.0 


4 


5 


2.0 


5 


10 


0.0 


6 


10 


0.5 


7 


10 


1.0 


8 


10 


2.0 


9 


10 


00 


10 


20 


0.5 


11 


20 


1.0 


12 


20 


2.0 


13 


00 


1.0 



Table 3: Abscissa key for Figures 5-16 

The corresponding plots for the errors in the analyzed wind 
fields are shown in Figures 8-10. The errors shown are the rms 
of the two winds, the values generally being quite close 
together, as can be seen by referring to Tables l.g.m in the 
Appendix 2. The errors in the analyzed wind fields generally 
follow the same trend as the errors in the analyzed height 
fields, except that the wind errors are more sensitive to the 
wind observation errors, this being especially prominent when 
there are no wind observations (at abscissa 9). One interesting 
thing to notice is that the curves for correlation functions 2-4 
are very close together, indicating an even greater dependence on 
the decay rate parameter, a, than the analyzed height errors 
exhibited. 

In order to assess the importance of a single observation 
and the effects of many missing observations, some 01 analysis 
errors were computed based on missing observed values. The 
missing observations are shown in Figures 1-3 for the three 
grids. The expected 01 height analysis errors are shown in 
Figures 11-13 for the three grids, and the expected 01 wind 
analysis errors in Figures 14-16. I will discuss Figure 11 in 
some detail, and a similar analysis follows for the other 
figures. 

The figure is slightly busy. The open circles and open 
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squares denote the values for correlation functions #4 and #1, 
respectively, the same values as given in Figure 5. The same 
symbol with the + denotes the analyzed values under the same 
assumptions but with observed values missing (appearing twice; 
once with observation of winds missing at the location noted 
previously, the second with observations of both winds and height 
missing). It can be noted that the missing wind observation is 
nearly undetectable while there is only a slight degradation with 
the height observation also missing. This is, however, more 
prominent in the case of smaller correlation in the upper graph 
for correlation function #1. The shaded symbols denote the 
values obtained when one-half (18) of the observations are 
missing (again, twice, once for missing winds, and once for 
missing heights and winds). At abscissa 9 the symbol overlays 
the others for which there are no wind observations and height 
observations at all locations. The nominal case (abscissa 7) 
shows about a 15 percent increase in the analysis errors when 18 
wind observations are missing, and about a 50 percent increase in 
the analysis errors for each correlation function when the entire 
observation is missing at 18 points. In this case, a 15 percent 
increase in the analysis errors corresponds to about a 6 percent 
decrease in "skill", where "skill" is taken as ( 1 - expected 
analysis error relative to background error), and the 50 percent 
increase in error to about a 20 percent decrease in skill. The 
general pattern of error for the various parameter values is 
generally the same as for that obtained for the entire 
observation set, the primary difference being for abscissa 9 
where the wind observations are all missing. 

The 01 analyzed height error for the EC grid (Figure 12) 
follows much the same pattern, except that the increase due to 
one missing observation is somewhat more significant, the total 
number of observations being 25 instead of 36 as in the MUS grid. 
Again, with one-half the observations missing, the error is 
increased by about 50 percent, and the drop in skill about 25 
percent . 

The MA grid results in Figure 13 differ since one missing 
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observation is nearly one-half of the total of three. Thus there 
are no shaded symbols in that figure, and a completely missing 
observation results in an increase of height analysis error of a 
few percent, but this again being about a 20 percent drop in 
skill . 

The plots for the wind analysis errors are given in Figures 
14-16. Note that some of symbols for separate correlation 
functions are overlaid in Figures 14 and 15. Of course, the wind 
analysis is sensitive to the loss of one or more observations, 
with the increase in wind analysis errors in Figure 14 for the 
MUS grid showing an increase of about 6-8 percent with 18 missing 
wind observations, and about 12-15 percent when the entire 
observation is missing at 18 points. The decrease in skill here 
corresponds to about 20 percent and 40 percent, respectively. On 
a relative basis the analysis errors are significantly larger for 
winds, with smaller relative increases in the error when 
observations are missing, however, in terms of skill level, the 
winds are more dependent on the observations. Again, the general 
character of the errors follows the same pattern, with the 
exception of abscissa 9, corresponding to no wind observations. 

The plots for the wind analysis errors for the EC grid in 
Figure 15 and the MA grid in Figure 16 reveal no surprises. The 
general pattern of Figure 15 is similar to Figure 14, while 
Figure 16 reveals that a very small skill is involved in this 
case, so missing observations have little affect. 

Unfortunately, none of the parameters varied above are 
really at the disposal of the practitioner. Still, the above 
information is a useful aid to understanding the 01 (and SI) 
process and how achievable accuracy is affected by the parameters 
in the process. 

The more important practical information is that given in 
Figures 17-22. Once again, the analysis errors are plotted 
versus a single abscissa which corresponds to various 
combinations of assumed parameter values. The nominal background 
rms error is 30 m for the height, as noted in the previous 
section, with the background wind errors depending on the 
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correlation function. Nominal height observation accuracy is 10 
m, and nominal wind observation accuracy is 1 m/sec. For a given 
spatial correlation function number, the assumed observation 
accuracies and spatial correlation function are varied. 
Observation accuracies vary by a factor of .5 or 2 from the 
nominal, while correlation function number varies by at most one. 

Consider Figure 17. The three "curves” for the MUS grid and 
spatial correlation functions 2, 3, and 4 will be discussed. As 
in previous figures, only the integer abscissa values have 
meaning, and the points are connected only to enable one to more 
easily see the effects of changing parameter values. It is 
immediately apparent that the parameters to which the SI scheme 
are most sensitive are embodied in abscissae 9, 12, 13, and to a 
lesser extent, 2, 4, and 5. Table 4 shows the relationship 
between the abscissae and the parameter variations, and we see 
that each of these abscissae except 9 are for assumed height 
observation errors that are twice the nominal value, and 
abscissae 9, 12, and 13 are for a misspecified correlation 
function (greater correlation). Abscissa -8, -9, -10, -11, -12, 
and -13 also show relatively larger SI errors, and each of these 
abscissa are for low assumed correlation as well as improper 
assumed height observation error. It appears that it is better 
to underestimate height observation errors than to overestimate 
them, although there is a peak (but smaller) at abscissa -2, 
where the height observation error is underestimated. 

Note that the graphs for correlation functions 3 and 4 are 
guite similar, while that for correlation function 2 differs 
somewhat for negative abscissae. Table 1 shows that the assumed 
correlation function (that is function 1) for large negative 
abscissae for correlation function 2 has a different decay rate, 
while for correlation functions 3 and 4, the decay rate is the 
same as that of 2 and 3. The relatively larger effect of the 
improper decay rate for the assumed correlation function is also 
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abscissa 


h error 


1 w error 


1 cor f . # 






nominal values 






10 


1 




m 






assumed values 




-13 


5 


0.5 




m-1 


-12 


5 


2 




m-1 


-11 


20 


0.5 




m-1 


-10 


20 


2 




m-1 


-9 


5 


1 




m-1 


-8 


20 


1 




m-1 


-7 


10 


0.5 




m-1 


-6 


10 


2 




m-1 


-5 


5 


0.5 




m 


-4 


5 


2 




m 


-3 


10 


1 




m-1 


-2 


5 


1 




m 


-1 


10 


0.5 




m 


0 


10 


1 




m 


1 


10 


2 




m 


2 


20 


1 




m 


3 


10 


1 




m+1 


4 


20 


2 




m 


5 


20 


0.5 




m 


6 


10 


2 




m+1 


7 


10 


0.5 




m+1 


8 


20 


1 




m+i 


9 


5 


1 




m+1 


10 


5 


0.5 




m+1 


11 


5 


2 




m+1 


12 


20 


0.5 




m+1 


13 


20 


2 




m+1 



Table 4: Abscissa key for Figures 17-22 

apparent in the SI error for correlation function 4 at abscissae 
9, 12, and 13, as noted above. Thus, the decay rate for the 
correlation function seems to be more important than the additive 
constant. 

For the EC grid, the results shown in Figure 18 indicate 
that the character of the three graphs is much the same. The 
outstanding difference is the significantly larger SI errors 
occurring for correlation function §2 when the assumed 
correlation function is §1 (abscissae -3 and -6 to -13). Again, 
the graph for correlation function #4 shows larger SI errors for 
assumed correlation function iS (larger abscissae). Both of 
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these cases correspond to misspecif ied decay rates. Looking 
higher to the MA grid, we see same effects: cases where the 

decay rate for the correlation function is misspecif ied yield 
larger increases in the SI height error than when the additive 
constant is misspecif ied. 

The results for the MA grid shown in Figure 19 imply that 
the most crucial parameter to have correct in such sparse regions 
is the decay rate for the correlation function, unfortunately the 
most difficult to estimate in such cases. Further, it appears it 
is probably best to underestimate the decay rate in sparse (or 
semi-sparse regions, such as EC) regions. In data dense regions 
the height errors are generally less sensitive to 
misspecif ication of the correlation function (with the exception 
of this being in combination with overestimates of the height 
observation error) . 

The corresponding plots for SI wind errors are in Figures 
20-22. Here the behavior of the errors seems to be less 
structured, with the smallest error often occurring for the 
correlation function corresponding to the least spatial 
correlation (SI height errors are generally a decreasing function 
of correlation function #). The behavior of the SI wind errors 
are also sensitive to the misspecif ication of the decay rate for 
the correlation function, as can be noted by the correlation 
function #2 values for large negative abscissae and for 
correlation function #4 for large positive abscissae. With that 
exception, the general behavior of the SI wind errors is much the 
same as for the SI height errors, with the primary dependence 
again being on the correct specification of the assumed height 
error. 

One additional bit of information can be squeezed from the 
data generated by this study. This concerns the relationship 
between the expected error based on the 01 calculations versus 
the actual expected error. Of course, since parameters are 
estimated, only the expected error for 01 can be calculated when 
expected errors are needed (e.g., see Goerss (1989)). We give 
three examples to show how this proceeds, and to show the 



14 



variation in the values. While the information can also be 
obtained by looking at the appropriate figures, the information 
is more precisely and as easily obtained from the tables in 
Appendix 1. 

( 1 ) Let us suppose that the assumed values of the 
observation errors are the twice the nominal ones and that the 
assumed spatial correlation function is #4. Height and wind 
observation errors equal 20 m and 2 m/sec, respectively. For 
the MUS grid the expected error for the heights (from Table 

1.1.4, or Figure 5, abscissa 12 on Cor Ftn 4) is 0.3462. If 
the actual values for the observation errors are 10 m for 
heights and 1 m for winds, with spatial correlation function 
#3, then the actual expected error for the MUS grid is (from 
Table 6. 1.3-4, or Figure 17, abscissa 13 on Cor Ftn 3) is 
0.3092, significantly smaller than the 01 calculation would 
indicate. 

(2) Now suppose the assumed values of the observation 
errors are one-half the nominal values and that the assumed 
correlation function is #3. Height and wind observation errors 
equal 5 m and 0.5 m/sec, respectively. For the MUS grid the 
expected error for the heights (from Table 1.1.3, or Figure 5, 
abscissa 2 on Cor Ftn 3) is 0.2291. Again if the correct 
values are the nominal values of 10 m and 1 m/sec for the 
height and wind observation errors, and the actual spatial 
correlation function is #4, then (from Table 6. 1.4-3, or Figure 
17, abscissa -13 on Cor Ftn 4 ) the correct expected error is 
0.2819. In this case the expected errors are significantly 
larger than the 01 calculation indicates. 

(3) In this case, suppose the spatial correlation function 
is correct, #4, and the assumed values of the observation 
errors are 5 m and 2 m/sec for heights and winds, respectively. 
Then the for the EC grid the expected height error (Table 

1.2.4, or Figure 6, abscissa 4 on Cor Ftn 4) is 0.3486. If the 
actual observation errors are 10 m and 1 m/sec for heights and 
winds, respectively, then the actual expected error (Table 

6. 2. 4-4, or Figure 18, abscissa -4 on Cor Ftn 4) is 0.3870; 
again the 01 expected error is smaller. 

As a general rule, expected error as calculated by 01 will 
be optimistic when the observation errors are underestimated or 
when the spatial correlation is overestimated. The latter 
usually has a greater influence. On the other hand, when 
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observation errors are overestimated or when the spatial 
correlation is underestimated, the expected error estimates 
computed by 01 will be pessimistic. Any expected error estimates 
from operational 01 sources should be treated with some caution, 
the examples above merely serving as an indication of the 
difficulties and not as a guide to the magnitude of the 
difference between the computed and actual values that may occur 
in practice. 

5 . 0 CONCLUDING REMARKS 

The results of this study demonstrate that SI analyzed 
height errors are more sensitive to the decay rate for the 
spatial correlation function than for the additive constant. 

While this study concentrated on the SOAR correlation function, 
similar results can be expected for other correlation functions 
which are controlled by parameters governing similar properties. 
The wind errors are even more sensitive to proper values for the 
decay rate, unsurprising since the wind correction is related to 
the derivative of the spatial correlation function for the 
heights . 

Another interesting observation is that it is better to 
"make the same mistake" relative to the observation error for the 
heights and winds. For example, if observed height error is 
specified as too large, smaller analysis errors occur when the 
observed wind errors are also specified as too large, rather than 
correct or too small. As a general rule, erring on the side of 
underestimating the observation error seems to result in smaller 
analysis errors than erring on the side of overestimation of the 
observation error. 

The effect of one missing wind observation is vanishingly 
small. This is true even when one observation constitutes a 
significant portion of the total amount of data, as in the MA 
grid case. However, two things come into play in this case to 
make the missing data still rather insignificant: (1) the 

missing data is close to another observation, and (2) the skill 
in this case is rather low anyway. Missing much data (about 
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one-half) shows significant decreases in skill, about 20-25 
percent for heights and up to 40 percent for the winds. 

The relationship between the expected error that can be 
computed using SI and the actual values of the expected error 
were explored briefly. Practitioners need to be cognizant of the 
fact that these two values may be significantly different from 
each other. 
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Figure 1: The Middle US (MUS) grid and observation locations. 

The +'s show the 2.5° grid. The open squares and circles and the 
shaded circle show observation locations. The shaded circle 
is a "missing observation" in one test run, and all circles are 
"missing observations" in another. 




Figure 2; The East Coast (EC) grid and observation locations. 
Symbols as in Figure 1. 
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Figxire 3: The Middle Atlantic (MA) grid and observation 

locations. Symbols as in Figure 1, except there is only a single 
"missing observation" . 
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Figure 4 : Spatial correlation functions. 
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Figure 5: Expected error in 01 height analysis under various 

observation error and spatial correlation function conditions for 
the MUS grid and observation set. 
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ERROR RELATIVE TO BACKGROUND ERROR 
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Figure 6; As in Figure 5, for the EC grid and observation set. 
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Figure 7: As in Figure 5, for the MA grid and observation set. 
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Figure 8: Expected error in 01 wind analysis under various 

observation error and spatial correlation function conditions for 
the MUS grid and observation set. 
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Figure 9: As in Figure 8, for the EC grid and observation set. 
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ERROR RELRTIVE TO BRCKGROUND ERROR 
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Figure 10: As in Figure 8, for the MA grid and observation set. 
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01 H ERROR N/ MISSING OBS, MUS 
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Figiire 11: Expected error in Ol height analysis under various 

observation error and spatial correlation function conditions for 
the MUS grid and observation set. The open square and circle 
mark the values for the entire observation set for correlation 
function numbers #1 and #4, respectively. The crossed square and 
circle mark two values; the lower for one missing wind 
observation, the upper one for both height and wind. The shaded 
square and circle mark two values; the lower one for 18 missing 
wind observations, the upper one for both heights and winds. 
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Figure 12: As in Figure 11, for EC grid and observation set, one 

and 12 missing observations. 
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Figure 13: As in Figure 11, for MA grid and observation set, one 

missing observation (only). 
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Figure 14: Expected error in 01 wind analysis under various 

observation error and spatial correlation function conditions for 
the MUS grid and observation set. Symbols as in Figure 11. 
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Figvire 15: As in Figure 14, for EC grid and observation set, one 

and 12 missing observations. 
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Figure 16: As in Figure 14, for MA grid and observation set, one 

missing observation (only). 
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HEIGHT SENSITIVITY PLOTS, MUS GRID 
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HEIGHT SENSITIVITY PLOTS, EC GRID 
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Figure 18: As in Figure 17 for the EC grid and observation 
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Figure 20: Expected error in the SI wind analysis under various 

assumed correlation function and parameter values (see Table 4 
for abscissae meanings), MUS grid and observation set. 
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Figure 22: As in Figure 20 for the MA grid and observation 



APPENDIX 1 : TABLES 



The tables giving the values plotted in the various figures, 
and some additional data as well, are given here. Most of the 
information is more readily accessible in the figures, however 
the tables are given for completeness. 

The numbering scheme for the tables is a key to the grid and 
observation set and the correlation function number used. The 
experiment whose results are given by the tables follows: 

Table l.g.m: Expected error under various values of parameters 

in 01, nominal observation data. 

g - the grid used, 1 for MUS, 2 for EC, 3 for MA. 
m - the spatial correlation function number, 1-5. 

Table 2.g.m: Expected error under various values of parameters 

in 01, wind observation missing at one point, as noted in text 
and Figures 1-3. g and m as in Table l.g.m, except m=l or 4. 

Table 3.g.m: Expected error under various values of parameters 

in 01, height and wind observation missing at one point, as noted 
in text and Figures 1-3. g and m as in Table l.g.m, except m=l 
or 4 . 

Table 4.g.in: Expected error under various values of parameters 

in 01, wind observations missing at about one-half the 
observation points, as noted in text and Figures 1-3. g and m as 
in Table l.g.m, except m=l or 4. 

Table 5.g.m: Expected error under various values of parameters 

in 01, height and wind observations missing at about one-half the 
observation points, as noted in text and Figures 1-3. g and m as 
in Table l.g.m, except m=l or 4. 

Table 6.g.m-n: Expected error under various values of parameters 

in SI. 

g - the grid used, 1 for MUS, 2 for EC, 3 for MA. 
m - the true spatial correlation function number, 1-5. 
n - the assumed spatial correlation function number, 1-5, 
but not differing more than one from m. 
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Obs . error 




Obs . 


error, heights 




winds 


0 

d 


5.0 


10.0 


20.0 


NOHolght 


0 

d 






0.4474 
O. 8018 
0.7978 






0.5 




0.4090 


0.4475 


0.5321 








0.7891 


O. 8019 


O. 8270 








0.7834 


0.7979 


0.8263 




0 

H 


0.3927 


0.4093 


0.4478 


0.5325 


0.7828 




0.7835 


0.7893 


0.8022 


O. 8273 


0.8739 




0.7771 


0.7837 


0.7982 


0.8266 


0.8814 


2.0 




0.4105 


0.4491 


0.5339 








O. 7905 


0.8034 


0.8285 








0.7847 


0.7992 


0.8277 




NoWlnds 






0.5434 
0.8911 
O. 8808 







Table 1.1.1: Expected error for observation accuracies 



Obs . error 
winds 


0.0 


Obs. 

5.0 


error, heights 

10.0 20.0 


NoHelght 


0.0 






0.2928 

0.6920 

0.6905 






0.5 




O. 2444 
0.6674 
0.6625 


0.2931 

0.6924 

0.6908 


O. 3742 
0.7240 
0.7277 




1.0 


O. 2174 
O. 6522 
0.6454 


0.2452 

0.6685 

0.6636 


0.2941 

0.6936 

0.6920 


0.3754 

0.7253 

0.7289 


0.6421 

0.7641 

0.7752 


2.0 




0.2483 

0.6728 

0.6677 


O. 2975 
0.6981 
0.6961 


O. 3800 
O. 7303 
0.7336 




NoMlndS 






O. 3794 
0.8108 
0.7970 







Table 1.1.2: Expected error for observation accuracies 



Obs . error 




Obs. 


error, heights 




winds 


O 

d 


5.0 


10.0 


20.0 


NoHelght 


0.0 






0.2775 

0.6957 

0.6947 






0.5 




0.2291 


0.2779 


0.3560 








0.6697 


0.6962 


0.7275 








0.6651 


0.6952 


0.7318 




0 

H 


O . 2007 


0.2300 


0.2789 


0.3574 


0.7078 




0.6525 


0.6710 


0.6976 


0.7291 


0.7645 




0.6457 


O. 6664 


0.6965 


0.7333 


0.7755 


;2.o 




0.2333 


0.2827 


0.3624 








0.6760 


0.7029 


0.7350 








0.6711 


0.7014 


0.7388 




NoWlnds 






0.3600 
0.8169 
O. 8032 







Table 1.1.3: Expected error for observation accuracies 
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Obs . error 




Obs . 


error , be i ghts 




vrinds 


0.0 


5.0 


10.0 


20.0 


NoHeight 


0-0 






0.2638 

0.6994 

0.6989 






0.5 




0.2157 


0.2642 


0.3391 








0.6721 


0.6999 


0.7307 








0.6678 


0.6994 


0.7355 




1.0 


0.1S59 


0.2166 


0.2653 


0.3406 


0.7571 




0.6529 


0.6737 


0.7016 


0.7326 


0.7649 




0.6460 


0.6692 


0.7010 


0.7373 


0,7759 


2.0 




0.2202 


0.2694 


0.3462 








0.6794 


0.7077 


0.7395 








0.6746 


0-7067 


0.7437 




NoWinds 






0.3427 

0.8230 

0.8095 







Table 1.1.4: Expect.ed error for observation accuracies 



Obs . error 
winds 


0.0 


Obs . 
5,0 


error, heights 

10.0 20.0 


NoHeight 


0.0 






O. 1983 
0.6255 
0.6283 






0.5 




0.1549 

0.5969 

0.5946 


0.1991 

0.6271 

0.6297 


0.2591 

0.6516 

0.6589 




1.0 


0.1183 

0.5653 

0.5583 


0.1567 

0.6008 

0.5982 


0.2015 

0.6315 

0.6338 


0.2625 

0.6567 

0.6636 


0.7283 

0.6772 

0,6879 


2.0 




0.1627 
0.6138 
O. 6101 


0.2095 

0.6462 

0.6474 


0.2740 

0.6740 

0.6798 




NoWinds 






0.2743 

0.7753 

0.7597 







Table 1.1.5: Expected error for observation accuracies 



Obs . error 




Obs . 


error, heights 




winds 


0 

d 


5.0 


10.0 


20.0 


NoHeight 


0 

6 






0.5903 

0.8543 

0.8543 






0.5 




0.5670 


0.5904 


0.6456 








0.8435 


0.8544 


0.8747 








0.8460 


0.8544 


0.8712 




0 

H 


0.5575 


0.5672 


0.5907 


0.6459 


0.8301 




0.8386 


0.8437 


0.8546 


0.8749 


0.9119 




0.8425 


0.8462 


0.8546 


0.8714 


,0 . 9058 


2.0 




0.5681 


0.5916 


0.6471 








0.8445 


0.8554 


0.8758 








0.8470 


0.8554 


0.8723 




Nowinds 






0.6631 

0.9121 

0.9190 







Table 1.2.1: Expected error for observation accuracies 
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Obs . error 
winds 


0.0 


Obs . 
5.0 


error, heights 

10.0 20.0 


NoHelght 


0.0 






0.4266 

0.7702 

0.7671 






0.5 




0.3972 

0.7515 

0.7516 


0.4269 

0.7705 

0.7674 


0.4854 

0.7939 

0.7886 




1.0 


0.3828 

0.7392 

0.7426 


0.3979 

0.7523 

0.7525 


0.4277 

0.7713 

0.7683 


0.4865 

0.7949 

0.7897 


0.7146 

0.8255 

0.8202 


2.0 




O . 4006 
0.7554 
0.7558 


0.4310 

0.7746 

0.7718 


0.4907 

0.7987 

0.7937 




NoWlnds 






0.5132 

0.8544 

0.8593 







Table 1.2.2: Expected error for obeervation accuracies 



Obs . error 
winds 


0.0 


Obs . 
5.0 


error, heights 

10.0 20.0 


NoHelght 


0.0 






0.3998 

0.7733 

0.7699 






0.5 




0.3696 

0.7537 

0.7535 


0.4001 

0.7737 

0.7703 


0.4588 

0.7969 

0.7915 




1.0 


0.3543 

0.7399 

0.7433 


0.3703 

0.7546 

0.7546 


0.4011 

0.7747 

0.7714 


0.4600 

0.7980 

0.7928 


0.7646 

0.8258 

0.8205 


2.0 




0.3733 

0.7582 

0.7584 


0.4046 

0.7786 

0.7756 


0.4647 

0.8025 

0.7975 




Nowinds 






0.4824 
0.8595 
O. 8640 







Table 1.2.3: Expected error for observation accuracies 



Obs . error 




Obs . 


error, heights 




winds 


0.0 


5.0 


10.0 


20.0 


NoHelght 


0.0 






0.3753 

0.7761 

0.7724 






0.5 




0.3445 


0.3756 


0.4336 








0.7556 


0.7765 


0.7994 








0.7552 


0.7729 


0,7940 




0 

H 


0.3284 


0.3454 


0.3767 


0.4350 


0.8030 




0.7402 


0.7568 


0.7778 


O . 8008 


0.8262 




0.7437 


0.7564 


0.7742 


0.7955 


0.8209 


2.0 




0.3486 


0.3806 


0.4402 








0.7609 


0.7823 


0.8061 








0.7609 


0.7790 


0.8010 




NoWlnds 






0.4544 

0.8643 

0.8684 







Table 1.2.4: Expected error for observation accuracies 
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obs . error 
vinds 


0.0 


Obs. 

5.0 


error, heights 

10.0 20.0 


NoHelght 


0.0 






0.2811 

0.7072 

0.7037 






0.5 




0.2498 

0.6865 

0.6853 


0.2819 

0.7084 

0.7050 


0.3378 

0.7277 

0.7236 




1.0 


0.2300 

0.6627 

0.6669 


0.2517 

0.6895 

0.6884 


0.2843 

0.7118 

0.7086 


0.3410 

0.7316 

0.7277 


0.7728 

0.7524 

0.7489 


2.0 




0.2580 

0.6995 

0.6990 


0.2925 

0.7234 

0.7207 


0.3519 

0.7453 

0.7417 




NoWlnds 






0.3686 

0.8249 

0.8274 







Table 1.2.5: Expected error for observation accuracies 



Obs . error 
vinds 


0.0 


Obs . 
5.0 


error, heights 

10.0 20.0 


NoHelght 


0.0 






0.9355 

0.9817 

0.9843 






0.5 




0.9320 

0.9806 

0.9837 


0.9355 

0.9817 

0.9843 


0.9454 

0.9842 

0.9859 




1.0 


0.9307 

0.9802 

0.9835 


0.9321 

0.9806 

0.9837 


0.9356 

0.9817 

0.9843 


0.9455 

0.9842 

0.9860 


0.9807 

0.9903 

0.9914 


2.0 




0.9323 

0.9808 

0.9838 


0.9358 

0.9818 

0.9844 


0.9457 

0.9843 

0.9861 




NoWlnds 






0.9520 

0.9896 

0.9922 







Table 1.3.1: Expected error for observation accuracies 



Obs . error 




Obs . 


error, heights 




vinds 


0.0 


5.0 


10.0 


20.0 


NoHelght 


0.0 






0.8567 

0.9585 

0.9649 






0.5 




0.8491 


0.8569 


0.8783 








0.9561 


0.9585 


0.9631 








0.9637 


0.9649 


0.9678 




O 

H 


0.8463 


0.8494 


0.8572 


0.8787 


0.9595 




0.9550 


0.9562 


0.9587 


0.9633 


0.9752 




0.9634 


0.9639 


0.9651 


0.9680 


0.9788 


2.0 




0.8508 


0.8586 


0.8801 








0.9569 


0.9594 


0.9641 








0.9646 


0.9658 


0.9687 




NoWlnds 






0.8934 

0.9767 

0.9846 







Table 1.3.2: Expected error for observation accuracies 
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Obs . error 
winds 


0.0 


Obs . 
5.0 


error, heights 

10.0 20.0 


NoHelght 


0.0 






0.8196 

0.9597 

0.9668 






0.5 




O. 8094 
0.9572 
0.9657 


0.8197 

0.9598 

0.9668 


0.8486 

0.9643 

0.9694 




1.0 


0.8058 

0.9561 

0.9655 


0.8098 

0.9574 

0.9659 


0.8201 
O . 9600 
0.9670 


0.8490 

0.9645 

0.9697 


0.9657 

0.9753 

0.9788 


2.0 




0.8112 

0.9582 

0.9667 


0.8216 

0.9609 

0.9678 


0.8505 

0.9654 

0.9705 




NoMlnds 






0.8533 

0.9783 

0.9866 







Table 1.3.3: Expected error for observation accuracies 



Obs. error 
winds 


0.0 


Obs . 
5.0 


error, heights 

10.0 20.0 


NoHelght 


0.0 






0.7769 

0.9607 

0.9680 






0.5 




0.7643 

0.9580 

0.9670 


0.7770 

0.9608 

0.9681 


0.8133 

0.9652 

0.9706 




1.0 


0.7598 

0.9567 

0.9667 


0.7647 

0.9583 

0.9672 


0.7775 

0.9610 

0.9683 


0.8137 

0.9655 

0.9709 


0.9708 

0.9753 

0.9789 


2.0 




0.7662 

0.9592 

0.9681 


0.7790 

0.9620 

0.9693 


0.8152 

0.9666 

0.9718 




NoWlnds 






0.8084 

0.9795 

0.9880 







Table 1.3.4: Expected error for observation accuracies 



Obs . error 
winds 


0.0 


Obs . 
5.0 


error, heights 

10.0 20.0 


NoHelght 


0.0 






0.6806 

0.9293 

0.9457 






0.5 




0.6626 

0.9254 

0.9444 


0.6811 

0.9296 

0.9460 


0.7332 

0.9373 

0.9497 




1.0 


0.6565 

0.9233 

0.9444 


0.6639 

0.9263 

0.9452 


0.6824 

0.9306 

0.9468 


0.7345 

0.9383 

0.9506 


0.9552 
0.9573 
, 0.9632 


2.0 




0.6684 

0.9294 

0.9482 


0.6870 

0.9341 

0.9499 


0.7389 

0.9420 

0.9537 




NoWlnds 






0.7311 

0.9642 

0.9818 







Table 1.3.5: Expected error for observation accuracies 
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Obs . error 




Obs . 


error, heights 




winds 


0.0 


5.0 


10.0 


20.0 


NoHelght 


0.0 






0.4493 

0.8064 

0.8005 






0.5 




0.4107 


0.4494 


0.5345 








0.7935 


0.8065 


0.8317 








0.7861 


0.8006 


0.8290 




1.0 


0.3944 


0.4110 


0.4497 


0.5349 


0.7879 




0.7879 


0.7938 


0.8067 


0.8320 


0.8794 




0.7798 


0.7864 


0.8008 


0.8293 


O. 8844 


2.0 




0.4123 


0.4510 


0.5363 








0.7949 


0.8078 


0.8331 








0.7874 


0.8019 


0.8303 




NoWlnds 






0.5434 

0.8911 

0.8808 







Table 2.1.1: Expected error for observation accxiracles 



Obs . error 




Obs . 


error , he Ights 




winds 


0.0 


5.0 


10.0 


20.0 


NoHelght 


0.0 






0.2652 

0.7040 

0.7018 






0.5 




0.2167 


0.2656 


0.3414 








0.6764 


0.7045 


0.7359 








O . 6706 


0.7023 


0.7386 




1.0 


0.1868 


0.2177 


0.2667 


0.3429 


0.7604 




0.6570 


0.6779 


0.7061 


0.7377 


0.7711 




0.6487 


0.6720 


0.7038 


0.7403 


0.7797 


2.0 




0.2212 


0.2708 


0.3484 








0.6834 


0.7120 


0.7444 








0.6772 


0.7094 


0.7467 




NoWlnds 






0.3427 

0.8230 

0.8095 







Table 2.1.4: Expected error for observation accuracies 



Obs . error 




Obs . 


error, heights 




winds 


0 

d 


5.0 


10.0 


20.0 


NoHelght 


0 

d 






0.5937 

0.8556 

0.8571 






0.5 




0 . 5702 


0.5938 


0.6497 








0.8448 


0.8557 


0.8761 








0.8487 


0.8571 


0.8742 




0 

H 


0.5606 


O . 5704 


0.5941 


0.6500 


0.8372 




O. 8399 


0.8450 


0.8559 


0.8763 


0.9138 




0.8451 


0.8489 


0.8574 


0.8745 • 


0.9096 


2.0 




0.5713 


0.5950 


0.6511 








0.8458 


0.8567 


0.8771 








0.8497 


0.8582 


O. 8753 




NoWlnds 






0.6631 

0.9121 

0.9190 







Table 2.2.1: Expected error for observation accuracies 
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Obs . error 
winds 


0.0 


Obs. 

5.0 


error, heights 

10.0 20.0 


NoHelght 


0.0 






0.3788 

0.7781 

0.7762 






0.5 




0.3474 

0.7574 

0.7585 


0.3792 

0.7785 

0.7766 


0.4384 

0.8019 

0.7984 




1.0 


0.3308 

0.7419 

0.7467 


0.3482 

0.7585 

0.7596 


0.3802 

0.7797 

0.7779 


0.4397 

0.8033 

0.7998 


0.8087 

0.8297 

0.8264 


2.0 




0.3513 

0.7626 

0.7640 


0.3840 

0.7842 

0.7826 


0.4448 

0.8085 

0.8052 




NoWinds 






0.4544 

0.8643 

0.8684 







Table 2.2.4: Expected error for observation accuracies 



Obs . error 
winds 


0.0 


Obs . 
5.0 


error, heights 

10.0 20.0 


NoHelght 


0.0 






0.9420 

0.9864 

0.9870 






0.5 




0.9382 

0.9854 

0.9863 


0.9421 

0.9865 

0.9870 


0.9525 

0.9891 

0.9888 




1.0 


0.9368 

0.9850 

0.9861 


0.9383 

0.9854 

0.9863 


0.9421 

0.9865 

0.9870 


0.9525 

0.9891 

0.9888 


0.9885 

0.9957 

0.9945 


2.0 




0.9384 

0.9854 

0.9864 


0.9423 

0.9865 

0.9871 


0.9527 

0.9891 

0.9889 




NoMinds 






0.9520 

0.9896 

0.9922 







Table 2.3.1: Expected error for observation accuracies 



Obs . error 




Obs . 


error, heights 




winds 


0.0 


5.0 


10.0 


20.0 


NoHelght 


O 

O 






0.8143 

0.9752 

0.9772 






0.5 




0.8014 


0.8144 


0.8515 








0.9741 


0.9753 


0.9783 








0.9765 


0.9773 


0.9795 




1.0 


0.7968 


0.8017 


O. 8148 


0.8518 


0.9805 




0.9738 


0.9743 


0.9754 


0.9785 


. 0.9869 




0.9764 


0.9767 


0.9775 


0.9797 


0.9862 


2.0 




0.8029 


0.8160 


0.8530 








0.9749 


0.9761 


0.9792 








0.9774 


0.9782 


0.9804 




NoMinds 






0.8378 

0.9887 

0.9914 







Table 2.3.4: Expected error for observation accuracies 
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Obs . error 
winds 


0.0 


Obs . 
5.0 


error, heights 

10.0 20.0 


NoHeight 


0.0 






0.4570 

0.8065 

0.8043 






0.5 




0.4198 

0.7941 

0.7904 


0.4571 

0.8066 

0.8044 


0.5400 

0.8310 

0.8320 




1 .0 


0.4042 

0.7888 

0.7843 


0.4202 

0.7944 

0.7906 


0.4574 

0.8069 

0.8047 


0.5404 

0.8312 

0.8323 


0,7892 

0.8772 

0.8864 


2.0 




0.4214 

0.7955 

0.7917 


0.4587 
O. 8080 
0.8057 


0.5418 

0.8324 

0.8333 




NoWinds 






0.5534 

0.8946 

0.8842 







Table 3.1.1: Expected error for observa-tion acctiracies 



Obs . error 
winds 


0.0 


Obs . 
5.0 


error, heights 

10.0 20.0 


NoHeight 


0.0 






0.2681 

0.7038 

0.7047 






0.5 




0.2207 

0.6771 

0.6744 


0.2685 

0.7043 

0.7052 


0.3439 

0.7348 

0.7411 




1.0 


0.1919 

0.6585 

0.6533 


0.2217 

0.6787 

0.6758 


0.2697 

0.7059 

0.7067 


0.3454 

0.7366 

0.7428 


0.7616 

0.7698 

0.7819 


2.0 




0.2253 

0.6844 

0.6811 


0.2738 

0.7120 

0.7123 


0.3509 

0.7435 

0.7490 




NoWinds 






0.3482 

0.8266 

0.8133 







Table 3.1.4: Expected error for observation accuracies 



Obs . error 
winds 


0.0 


Obs . 
5.0 


error, heights 

10.0 20.0 


NoHeight 


0.0 






O . 6000 
0.8595 
0.8592 






0.5 




0.5771 
O. 8496 
0.8512 


0.6001 

0.8596 

0.8593 


0.6544 

0.8784 

0.8756 




1.0 


0.5678 

0.8452 

0.8478 


0.5773 
0.8498 
O. 8514 


O . 6003 
0.8598 
0.8595 


0.6547 

0.8786 

0.8758 


0.8372 

0.9138 

0.9096 


2.0 




0.5782 

0.8506 

0.8522 


0.6013 

0.8606 

0.8603 


0.6558 

0.8794 

0.8767 




NoWinds 






0.6723 

0.9166 

0.9221 







Table 3.2.1: Expected error for observation accuracies 
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Obs. error 
winds 


0.0 


Obs. 

5.0 


error, heights 

10.0 20.0 


NoHeiqht 


0.0 






0.3816 

0.7805 

0.7778 






0.5 




0.3509 

0.7614 

0.7609 


0.3820 

0.7809 

0.7782 


0.4406 

0.8029 

0.7991 




1.0 


0-3348 

0.7474 

0.7498 


0.3518 

0.7625 

0.7621 


0.3831 

0.7821 

0.7795 


0.4420 

0.8043 

0.8006 


0.8087 

0.8297 

0.8264 


2.0 




0.3550 

0.7667 

0.7665 


0.3870 

0.7867 

0.7843 


0.4472 

0.8096 

0.8060 




NoWinds 






0.4611 

0.8689 

0.8721 







Table 3.2.4: Expected error for observation accuracies 



Obs . error 
winds 


0.0 


Obs. 

5.0 


error, heights 

10.0 20.0 


NoHeight 


0.0 






0.9569 

0.9902 

0.9896 






0.5 




0.9543 

0.9897 

0.9892 


0.9569 

0.9902 

0.9896 


0.9643 

0.9914 

0.9908 




1.0 


0.9534 

0.9896 

0.9891 


0.9544 

0.9897 

0.9893 


0.9570 

0.9902 

0.9896 


0.9643 

0.9915 

0.9908 


0.9885 

0.9957 

0.9945 


2.0 




0.9545 

0.9898 

0.9893 


0.9572 

0.9902 

0.9897 


0.9645 

0.9915 

0.9909 




NoWinds 






0.9688 

0.9945 

0.9951 







Table 3.3.1: Expected error for observation accuracies 



Obs . error 




Obs . 


error, heights 




winds 


0.0 


5.0 


10.0 


20.0 


NoHeight 


O 

d 






0.7889 

0.9695 

0.9747 






0.5 




0.7747 


0.7890 


0.8264 








0.9659 


0.9695 


0.9755 








0.9732 


0.9747 


0.9777 




O 

H 


0.7691 


0.7750 


0.7893 


0.8267 


0.9805 




0.9642 


0.9661 


0.9697 


0.9757 


0.9869 




0.9727 


0.9734 


0.9749 


0.9779 


0.9862 


i.o 




0.7761 


0.7903 


0.8277 








0.9666 


0.9702 


0.9762 








O . 9741 


0.9756 


0.9786 




NoWinds 






0.8084 

0.9795 

0.9880 







Table 3.3.4: Expected error for observation accuracies 
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Obs. error 




Obs . 


error, heights 




winds 


0.0 


5.0 


10.0 


20.0 


NoHeight 


0.0 






0.4933 

0.8457 

0.8406 






0.5 




0.4540 


0.4933 


0.5838 








O. 8321 


0.8457 


0.8734 








0.8253 


0.8407 


0.8720 




1.0 


0.4377 


0.4542 


0.4935 


0.5840 


0,8774 




0.8262 


O. 8322 


0.8459 


0.8735 


0.9305 




0.8187 


0.8255 


0.8408 


0.8721 


0.9399 


2.0 




0.4548 


0.4942 


0.5848 








0.8328 


0.8465 


0.8741 








0.8260 


0.8413 


0.8727 




NoMinds 






0.5434 

0.8911 

0.8808 







Table 4.1.1: Expected error for observation accuracies 



Obs. error 
winds 


O 

d 


Obs . 
5.0 


error, heights 

10.0 20.0 


NoHeight 


O 

d 






0.3014 

0.7570 

0.7534 






0.5 




0.2477 

0.7244 

0.7167 


0.3016 

0.7573 

0.7537 


0.3920 

0.7983 

0.8016 




1.0 


0.2167 

0.7026 

0.6924 


0.2482 

0.7253 

0.7175 


0.3022 

0.7583 

0.7545 


0.3929 

0.7994 

0.8025 


0.8412 

0.8541 

0.8688 


2.0 




0.2500 

0.7285 

0.7202 


0.3045 

0.7617 

0.7574 


0.3962 
0.8035 
O . 8060 




NoWinds 






0.3427 

0.8230 

0.8095 







Table 4.1.4: Expected error for observation accuracies 



Obs . error 




Obs - 


error, heights 




winds 


o 

d 


5.0 


10.0 


20.0 


NoHeight 


o 

d 






0.6297 

0.8847 

0.8860 






0.5 




0.6038 


0.6297 


0.6915 








0.8731 


0.8847 


0.9076 








0.8768 


0.8860 


0.9053 




O 

H 


0.5932 


0.6039 


0.6298 


0.6916 


0.9039 




0.8679 


0.8732 


0.8848 


0.9077 


0.9537 




0.8729 


0.8769 


0.8862 


0.9054 


0.9481 


2.0 




0.6043 


0.6303 


0.6922 


' 






0.8735 


0.8852 


0.9081 








0.8774 


0.8866 


0.9059 




NoWinds 






0.6632 

0.9121 

0.9190 







Table 4.2.1: Expected error for observation accuracies 
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Obs • error 




Obs. 


error, heights 




winds 


0.0 


5.0 


10.0 


20.0 


NoHeight 


0.0 






0.4211 
O. 8237 
O. 8191 






O. 5 




0.3826 


0.4212 


0.4911 








0.7977 


0.8239 


0.8555 








0.7970 


0.8193 


O. 8483 




1 .O 


0.3617 


0.3830 


0.4217 


0.4918 


0.8703 




0.7796 


0.7982 


0.8245 


0.8562 


0.8978 




0.7826 


0.7976 


0.8200 


0.8491 


0.8897 


2.0 




0.3843 


0.4234 


0.4945 








O. 8003 


0.8268 


0.8590 








O. 7999 


0.8226 


O. 8522 




NoWinds 






0.4544 

0.8643 

0.8684 







Table 4.2.4: Expect.ed error for obaervat.lon accuracies 
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Obs . error 




Obs . 


error t heights 




vinds 


0.0 


5.0 


10.0 


20.0 


NoHeight 


0 

d 






0.6457 

0.8893 

0.8900 






0.5 




0.6247 


0.6458 


0.6989 








0.8831 


0.8894 


0.9024 








0.8838 


0.8901 


0.9042 




1.0 


0.6166 


0.6250 


0.6461 


0.6991 


O. 8774 




0.8806 


0.8833 


0.8895 


0.9026 


0.9305 




O. 8814 


0.8840 


0.8902 


0.9043 


0.9399 


2.0 




0.6258 


0.6470 


0.7001 








0.8839 


0.8902 


0.9033 








0.8845 


0.8908 


0.9049 




NoMinds 






0.7214 

0.9415 

0.9390 







Table 5.1.1: Expected error for observation accuracies 



Obs . error 
winds 


0.0 


Obs . 
5.0 


error, heights 

10.0 20.0 


NoHeight 


0.0 






O . 4067 
0.8070 
0.8035 






0.5 




0.3747 

0.7926 

0.7875 


0.4070 

0.8074 

0.8039 


0.4729 
O. 8263 
O. 8289 




1.0 


0.3595 

0.7838 

0.7789 


0.3756 

0.7936 

0.7885 


0.4080 

0.8085 

0.8049 


0.4741 

0.8275 

0.8300 


0.8412 
0.8541 
O. 8688 


2.0 




0.3789 

0.7975 

0.7920 


0.4118 

0.8126 

0.8085 


0.4787 

0.8322 

0.8340 




NoWinds 






0.4887 

0.8949 

0.8862 







Table 5.1.4: Expected error for observation accuracies 



Obs . error 
winds 


0.0 


Obs . 
5.0 


error, heights 

10.0 20.0 


NoHeight 


0,0 






0.7397 
O. 9175 
0.9164 






0.5 




0.7260 

0.9123 

0.9121 


0.7398 

0.9175 

0.9164 


0.7753 

0.9286 

0.9260 




1.0 


0.7207 

0.9103 

0.9103 


0.7262 

0.9125 

0.9122 


0.7400 

0.9177 

0.9166 


0.7755 

0.9287 

0.9261 


0.9039 

0.9537 

0.9481 


2.0 




0.7268 

0.9129 

0.9127 


0.7406 

0.9181 

0.9171 


0.7762 

0.9292 

0.9266 




NoWinds 






0.7946 

0.9523 

0.9558 







Table 5.2.1: Expected error for observation accuracies 
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Obs . error 




Obs . 


error, heights 




winds 


0.0 


5.0 


10.0 


20.0 


NoHelght 


0.0 






0.5221 

0.8580 

0.8540 






0.5 




0.5022 


0.5224 


0.5674 








0.8465 


0.8582 


0.8738 








0.8440 


0.8543 


0.8682 




1.0 


0.4931 


0.5028 


0.5232 


0.5684 


0.8703 




0.8396 


0.8472 


0.8590 


0.8746 


0.8978 




0.8382 


0.8448 


0.8552 


0.8691 


0.8897 


2.0 




0.5054 


0.5262 


0.5722 








0.8498 


0.8618 


0.8779 








0.8477 


0.8583 


0.8727 




NoWinds 






0.5899 

0.9181 

0.9214 







Table 5.2.4: Expected error for obaervatlon accuracies 
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Assmd Obs. 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.4531 

0.8052 

0.8015 


0.4478 

0.8022 

0.7982 


0.4717 

0.8118 

0.8087 


1.0 


0.4531 

0.8052 

0.8015 


0.4478 

0.8022 

0.7982 


0.4716 
O. 8117 
0.8087 


2.0 


0.4532 

0.8052 

0.8015 


0.4478 

0.8022 

0.7982 


0.4713 

0.8117 

0.8087 



Table 6. 1.1—1: Variation of expected error 



Assmd 


Obs. 


Assumed 


Obs . error , 


heights 


err, winds 


5.0 


10.0 


20.0 


0.5 




0.4583 


0.4709 


0.5476 






0.8064 


0.8120 


0.8379 






0.8020 


0.8081 


0.8371 


1.0 




0.4580 


0.4699 


0.5457 






0.8061 


0.8116 


O. 8374 






0.8017 


0.8077 


0.8366 


2.0 




0.4568 


0.4666 


0.5387 






0.8054 


0.8105 


0.8356 






0.8011 


0.8067 


O. 8349 


Table 6 


> . 1 . 1> 


-2: Variation of expected error 


Assmd 


Obs. 


Assumed 


Obs. error , 


heights 


err, winds 


5.0 


10.0 


20.0 


0.5 




0.3200 


0.3096 


0.3242 






0.7226 


0.7068 


0.7032 






0.7232 


0.7054 


0.7014 


1.0 




0.3201 


0.3097 


0.3243 






0.7227 


0.7069 


0.7033 






0.7233 


0.7055 


0.7015 


2.0 




0.3205 


0.3101 


0.3247 






0.7230 


0.7073 


0.7036 






0.7235 


0.7058 


0.7018 


Table 6 


1.1.2- 


-l: Variation of expected error 


Assmd 


Obs . 


Assumed 


Obs. error. 


heights 


err, winds 


5.0 


10.0 


20.0 


0.5 




0.3064 


0.2941 


0.3217 






0.7060 


0.6936 


0.7094 






0.7054 


0.6920 


0.7100 


o 

H 




0.3066 


0.2941 


0.3212 






0.7060 


, 0.6936 


0.7093 






0.7055 


0.6920 


0.7099 


2.0 




0.3076 


0.2942 


0.3197 






0.7066 


0.6938 


0.7090 






0.7060 


0.6921 


0.7096 



Table 6. 1.2-2: Variation of expected error 



Assmd Obs. 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.3045 


0.2945 


0.3279 




0.7038 


0.6939 


0.7127 




0.7030 


0.6923 


0.7139 


0 

H 


0.3047 


0.2944 


0.3273 




0.7039 


0.6939 


0.7126 




0.7031 


0.6922 


0.7137 


2.0 


0.3059 


0.2944 


0.3255 




0.7045 


0.6941 


0.7123 




0.7037 


0.6924 


0.7134 



Table 6. 1.2-3: Variation of expected error 



Assmd 


Obs . 


Assumed 


Obs . error , 


heights 


err, winds 


5.0 


10.0 


20.0 


0.5 




0.2945 


0.2793 


0.3016 






0.7149 


0.6980 


0.7108 






0.7152 


0.6968 


0.7116 


O 

H 




0.2947 


0.2793 


0.3012 






0.7150 


0.6980 


0.7107 






0.7153 


0.6968 


0.7115 


2.0 




0.2958 


0.2796 


0.2999 






0.7156 


0.6982 


0.7104 






0.7159 


0.6971 


0.7112 


Table 6 


.1.3- 


-2: Variation of expected error 


Assmd 


Obs. 


Assumed 


Obs . error , 


heights 


err, winds 


5.0 


10.0 


20.0 


0.5 




0.2922 


0.2789 


0.3059 






0.7121 


0.6977 


0.7139 






0.7122 


0.6965 


0.7152 


1.0 




0.2925 


0.2789 


0.3054 






0.7122 


0.6976 


0.7137 






0.7123 


0.6965 


0.7150 


2.0 




0.2937 


0.2791 


0.3038 






0.7129 


0.6979 


0.7134 






0.7130 


0.6967 


0.7147 


Table 6 


.1.3-3: Variation of expected error 


Assmd 


Obs . 


Assixmed 


Obs. error. 


heights 


err, winds 


5.0 


10.0 


20.0 


0.5 




0.2901 


0.2793 


0.3118 






0.7096 


0.6979 


0.7170 






0.7095 


0.6968 


0.7188 


1.0 




, O . 2904 


0.2792 


0.3112 






0.7097 


0.6979 


0.7168 






0.7096 


0.6968 


0.7186 


2.0 




0.2918 


0.2792 


0.3092 






0.7106 


0.6982 


0.7165 






0.7104 


0.6970 


0.7183 



Table 6.1. 3-4: Variation of expected error 
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Assmd Obe . 
err, winds 


Assumed 

5.0 


Obe . error , 
10.0 


heights 

20.0 


0.5 


0.2819 

0.7213 

0.7224 


0.2655 

0.7018 

0.7013 


0.2868 

0.7151 

0.7166 


1.0 


0.2B22 

0.7214 

0.7225 


0.2656 

0.7018 

0.7012 


0.2864 

0.7150 

0.7164 


2.0 


0.2835 

0.7222 

0.7232 


0.2659 

0.7021 

0.7015 


0.2850 

0.7147 

0.7161 



Tablo 6. 1.4-3: Variation of oxpoctod error 



Assmd Obs . 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.2796 

0.7182 

0.7190 


0.2653 

0.7016 

0.7010 


0.2917 

0.7181 

0.7200 


1.0 


0.2799 

0.7184 

0.7192 


0.2653 

0.7016 

0.7010 


0.2912 

0.7179 

0.7198 


2.0 


0.2814 

0.7193 

0.7201 


0.2655 

0^7019 

0.7013 


0.2894 

0.7176 

0.7195 


Table 6.1.4- 


-4: Variation of expected error 


Assmd Obs . 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.2693 

0.7051 

0.7046 


0.2777 

0.7105 

0.7107 


0.3320 

0.7369 

0.7420 


1.0 


0.2695 

0.7050 

0.7046 


0.2760 

0.7097 

0.7100 


0.3287 

0.7358 

0.7408 


2.0 


0.2716 

0.7072 

0.7067 


0.2724 

0.7092 

0.7095 


0.3190 

0.7341 

0.7391 


Table 6. 1.4-5: Variation of expected error 


Assmd Obs. 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.2383 

0.7005 

0.7084 


0.2095 

0.6431 

0.6465 


0.2117 

0.6370 

0.6398 


1.0 


0.2388 

0.7009 

0.7088 


0.2099 

0.6434 

0.6468 


0.2117 

0.6371 

0.6399 


2.0 


0.2408 

0.7029 

0.7106 


0.2113 

0.6448 

0.6481 


0.2120 

0.6377 

0.6404 



Table 6. 1.5-4: Variation of expected error 



Assmd Obs . 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.2187 

0.6587 

0.6640 


O . 2016 
0.6316 
0.6339 


0.2237 

0.6467 

0.6517 


1.0 


0.2198 

0.6595 

0.6647 


0.2015 

0.6315 

0.6338 


0.2224 

0.6461 

0.6511 


2.0 


0.2242 

0.6645 

0.6694 


0.2024 

0.6331 

0.6354 


0.2192 

0.6462 

0.6511 


Table 6.1. 5-! 


a: Variation of expected error 


Assmd Obs. 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.5939 

0.8574 

0.8564 


0.5907 

0.8546 

0.8546 


0.6057 

0.8626 

0.8608 


1.0 


0.5940 

0.8574 

0.8564 


0.5907 

0.8546 

0.8546 


0.6056 

0.8626 

0.8607 


2.0 


0.5941 

0.8574 

0.8565 


0.5907 

o,a&46 

0.8546 


0.6054 

0,8626 

0.8607 


Table 6.2.1- 


-1 : Variation of expected error 


Assmd Obs . 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.6311 

0.8631 

0.8627 


0.6368 

0.8681 

0.8663 


0.6817 

0.8884 

0.8832 


O 

H 


0.6304 

0.8628 

0.8625 


0.6356 

0.8677 

0.8659 


0.6797 

0.8879 

0.8827 


2.0 


0.6279 

0.8620 

0.8617 


0.6313 

0.8665 

0.8647 


0.6725 

0.8863 

0.8810 


Table 6.2. 1-2: Variation of expected error 


Assmd Obs. 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.4682 

0.8059 

0.7957 


0.4655 

0.7913 

0.7857 


0.4823 

0.7865 

0.7826 


1.0 


. 0.4683 

O. 8060 
0.7957 


0.4656 

0.7913 

0.7858 


0.4824 

0.7865 

0.7826 


2.0 


0.4689 
O . 8063 
0.7961 


0.4661 

0.7916 

0.7861 


0.4830 

0.7868 

0.7829 



Table 6. 2.2-1: Variation of expected error 
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Assmd Obs . 
err, winds 


Assumed 

5.0 


Obs . error , 
10.0 


heights 

20.0 


0.5 


0.4351 

0.7818 

0.7759 


0.4278 

0.7713 

0.7683 


0.4458 

0.7830 

0.7785 


1.0 


0.4353 

0.7818 

0.7760 


0.4277 

0.7713 

0.7683 


0.4455 

0.7829 

0.7784 


2.0 


0.4362 

0.7823 

0.7765 


0.4279 

0.7715 

0.7685 


0.4445 

0.7827 

0.7782 


Table 6.2.2- 


-2: Variation of expected error 


Assmd Obs . 


Assumed 


Obs . error , 


heights 


err, winds 


5.0 


10.0 


20.0 


0.5 


0.4348 

0.7800 

0.7748 


0.4287 

0.7717 

0.7687 


0.4499 

0.7856 

0.7809 


1.0 


0.4350 

0.7801 

0.7749 


0.4287 

0.7717 

0.7687 


0.4495 

0.7855 

0.7808 


2.0 


0.4360 

0.7807 

0.7754 


0.4287 

0.7719 

0.7689 


0.4483 

0.7853 

0.7806 



Tabid 6. 2. 2-3: Variation of oxpected error 



Assmd Obs . 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.4111 

0.7895 

0.7824 


0.4025 

0.7753 

0.7720 


0.4193 

0.7848 

0.7803 


0 

H 


0.4113 

0.7896 

0.7825 


O. 4025 
0,7753 
0.7720 


0.4191 

0.7847 

0.7802 


2.0 


0.4123 

0.7902 

0.7830 


0.4028 

0.7755 

0.7722 


0.4183 

0.7845 

0.7800 



Tabid 6.2. 3-2: Variation of expected error 



Assmd Obs. 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.4091 

0.7868 

0.7803 


0.4011 

0.7747 

0.7714 


0,4192 

0.7867 

0.7820 


1.0 


0.4093 

0.7869 

0.7804 


0.4011 

0.7747 

0.7714 


0.4188 

0.7866 

0.7818 


2.0 


0.4104 

0.7875 

0.7810 


0.4012 

0.7749 

0.7716 


0.4178 

0.7863 

0.7816 



Table 6.2. 3-3: Variation of expected error 



Assmd Obs . 
err, winds 


Assumed 

5.0 


Obs . error , 
10.0 


heights 

20.0 


0.5 


0.4081 

0.7846 

0.7788 


0.4015 

0.7750 

0.7717 


0.4230 

0.7890 

0.7841 


0 

H 


0.4084 

0.7847 

0.7789 


0.4014 

0.7749 

0.7716 


0.4226 

0.7889 

0.7840 


2.0 


0.4096 

0.7855 

0.7797 


0.4015 

0.7752 

0.7719 


0.4213 

0.7886 

0.7838 



Table 6. 2. 3-4: Variation of expected error 



Assmd Obs. 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.3868 

0.7942 

0.7864 


0.3771 

0.7780 

0.7745 


0.3924 

0.7879 

0.7832 


1.0 


0.3871 

0.7943 

0.7866 


0.3771 

0.7780 

0.7744 


0.3921 

0.7877 

0.7830 


2.0 


0.3883 

0.7951 

0.7873 


0.3774 

0.7783 

0.7747 


0.3912 

0.7875 

0.7828 


Table 6.2.4- 


-3: Variation of expected error 


Assmd Obs . 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.3854 

0.7914 

0.7844 


0.3767 

0.7778 

0.7742 


0.3948 

0.7899 

0.7851 


O 

H 


0.3857 

0.7916 

0.7846 


0.3767 

0.7778 

0.7742 


0.3944 

0.7898 

0.7849 


2.0 


0.3870 

0.7925 

0.7854 


0.3769 

0.7780 

0.7745 


0.3933 

0.7896 

0.7847 


Table 6. 2. 4-4: Variation of expected error 


Assmd Obs . 
err, winds 


Assumed 

5.0 


Obs . error , 
10.0 


heights 

20.0 


0.5 


0.3872 

0.7834 

0.7795 


0.3908 

0.7872 

0.7828 


0.4253 

0.8060 

0.8003 


1.0 


• 0.3869 

0.7831 
0.7793 


0.3891 

0.7864 

0.7820 


0.4223 

0.8049 

0.7993 


2.0 


0.3872 

0.7841 

0.7803 


0.3849 

0.7852 

0.7811 


0.4137 
O. 8029 
0.7975 



Table 6. 2. 4-5: Variation of expected error 
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Aesmd Obs . 
err, winda 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 




0-5 


O . 8889 


0.8913 


0.9013 






0.9678 


0.9669 


0.9670 






0.9740 


0.9732 


0.9723 




O 

H 


0.8890 


0.8914 


0.9014 






0.9678 


0.9669 


0.9671 






0,9740 


0.9732 


0.9724 




2.0 


O. 8892 


0.8916 


0.9016 






0.9679 


0.9670 


0.9672 






0.9741 


0-9733 


0.9724 





Table 6. 3. 2-1: Variation of expected error 



Assmd Obs. 
err, winds 


Assumed 

5.0 


Obs . error , 
10.0 


heights 

20.0 


0.5 


0.3106 

0.7703 

0.7549 


0.2946 

0.7240 

0.7196 


O. 3036 
0.7203 
0.7167 


1.0 


0.3112 

0.7708 

0.7553 


0.2950 

0.7243 

0.7199 


0.3038 

0.7204 

0.7168 


2.0 


0,3136 

0.7729 

0.7574 


0.2967 

0.7258 

0.7213 


O. 3049 
0.7212 
0-7176 



Table 6. 2.5-4: Variation of expected error 



Assmd Obs. 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.8581 

0.9595 

0.9653 


0.8572 

0.9587 

0.9651 


0.8620 

0.9603 

0.9658 


O 

H 


0.8581 

0.9595 

0.9653 


0.8572 

0.9587 

0.9651 


O. 8619 
0.9603 
0.9658 


2.0 


0.8583 

0.9596 

0.9653 


0.8573 

0.9587 

0.9651 


0.8620 

0.9603 

0.9658 



Table 6. 3. 2-2: Variation of expected error 



Assmd Obs . 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.2952 

0.7322 

0.7253 


0.2844 

0.7119 

0.7087 


0.3019 

0.7231 

0.7192 


1.0 


0.2961 

0,7329 

0.7260 


0.2843 

0.7118 

0.7086 


0.3009 

0.7226 

0.7188 


2.0 


0.3001 

0.7371 

0.7300 


0.2852 

0.7132 

0.7100 


0.2990 

0.7228 

0.7190 


Table 6.2.5- 


-5: Variation of expected error 


Assmd Obs . 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.9359 

0.9819 

0.9844 


0.9356 

0.9817 

0.9843 


0.9379 

0.9825 

0.9847 


1.0 


0.9359 

0.9819 

0.9844 


0.9356 

0.9817 

0.9843 


0.9379 

0.9825 

0.9847 


2.0 


0.9359 

0.9819 

0.9844 


0.9356 

0.9817 

0.9843 


0.9379 

0.9825 

0.9847 


Table 6.3. 1-1: Variation of expected error 


Assmd Obs . 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.9916 
O. 9878 
0.9901 


0.9873 

0.9879 

0.9902 


0.9807 

0.9897 

0.9908 


1.0 


0.9907 

0.9876 

0.9899 


0.9863 

0.9878 

0.9900 


0.9796 

0.9895 

0.9907 


2.0 


0.9872 

0.9871 

0.9894 


0.9826 

0.9872 

0.9895 


0.9757 

0.9889 

0.9901 



Assmd Obs . 
err, winds 


Assumed 

5.0 


Obs . error , 
10.0 


heights 

20.0 


0.5 


0.8680 

0.9595 

0.9656 


O. 8649 
0.9590 
0.9656 


0.8639 

0.9611 

0.9667 


1 .O 


0.8680 

0.9595 

0.9656 


0.8649 

0.9590 

0.9656 


0.8638 

0.9611 

0.9667 


2.0 


0.8681 

0.9596 

0.9657 


0.8649 

0.9590 

0.9656 


0.8639 

0.9611 

0.9668 



Table 6. 3. 2-3: Variation of expected error 



Assmd Obs . 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


O. 5 


0.8284 

0.9616 

0.9683 


0.8294 

0.9604 

0.9677 


0.8402 

0.9613 

0.9675 


1.0 


0.8284 
0.9616 
• 0.9683 


O. 8294 
O . 9604 
0.9677 


0.8402 

0.9613 

0.9675 


2.0 


0.8286 

0.9617 

0.9683 


0.8295 

0.9604 

0.9677 


0.8402 

0.9613 

0.9676 



Table 6. 3. 3-2: Variation of expected error 



Table 6. 3. 1-2: Variation of expected error 
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Asamd Obs . 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.8212 

0.9610 

0.9673 


0.8201 

0.9600 

0.9671 


0.8261 

0.9617 

0.9677 


1.0 


0.8212 

0.9610 

0.9673 


0.8201 

0.9600 

0.9670 


0.8261 

0.9617 

0.9677 


2 .O 


0.8213 

0.9611 

0.9673 


0.8202 

0.9601 

0.9671 


0.8261 

0.9617 

0.9678 



Table 6. 3. 3-3: Variation of expected error 



Assmd Obs . 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.8266 

0.9609 

0.9674 


0.8238 

0.9602 

0.9673 


0.8249 

0.9623 

0.9683 


1.0 


0.8266 

0.9609 

0.9674 


0.8238 

0.9602 

0.9673 


0.8249 

0.9623 

0.9683 


2.0 


0.8268 

0.9610 

0.9675 


0.8238 

0.9602 

0.9674 


0.8249 

0.9623 

0.9684 


Table 6.3.3- 


-4: Variation of expected error 


Assmd Obs . 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.7810 

0.9626 

0.9692 


0.7818 

0.9612 

0.9687 


0.7940 

0.9624 

0.9689 


O 

H 


0.7810 

0.9626 

0.9692 


0.7818 

0.9612 

0.9687 


0.7940 

0.9624 

0.9689 


2.0 


0.7812 

0.9627 

0.9692 


0.7819 

0.9613 

0.9687 


0.7940 

0.9624 

0.9689 


Table 6. 3. 4-3: Variation of expected error 


Assmd Obs . 
err, winds 


Assumed 

5.0 


Obs. error, 
10.0 


heights 

20.0 


0.5 


0.7787 

0.9621 

0.9686 


0.7775 

0.9610 

0.9684 


0.7845 

0.9628 

0.9691 


O 

H 


0.7787 

0.9622 

0.9686 


0.7775 

0.9610 

0.9683 


0.7845 

0.9627 

0.9691 


2.0 


0.7789 

0.9623 

0.9687 


0.7776 

0.9611 

0.9684 


0.7846 

0.9628 

0.9691 



Assmd 


Obs . 


Assumed 


Obs. error. 


heights 


err, winds 


5.0 


10.0 


20.0 


0.5 




0.8047 


0.8001 


0.7985 






0.9668 


0.9667 


0.9684 






0.9733 


0.9734 


0.9743 


o 

H 




0.8030 


0.7983 


0.7967 






0.9664 


0.9662 


0.9679 






0.9728 


0.9729 


0.9739 


2.0 




0.7984 


0.7930 


0.7912 






0.9653 


0.9647 


0.9663 






0.9717 


0.9717 


0.9727 


Table < 


6. 3.4-5: Variation of expected error 


Assmd 


Obs . 


Assumed 


Obs . error , 


heights 


err, winds 


5.0 


10.0 


20.0 


0.5 




0.6957 


0.6980 


0.7169 






0.9403 


0.9370 


0.9387 






0.9536 


0.9524 


0.9522 


o 

H 




0.6960 


0.6983 


0.7172 






0.9405 


0.9372 


0.9388 






0.9538 


0.9525 


0.9523 


2.0 




0.6972 


0.6995 


0.7184 






0.9412 


0.9379 


0.9394 






0.9544 


0.9531 


0.9529 


Table < 


6. 3. 5-4: Variation of expected error 


Assmd 


Obs . 


Assumed 


Obs. error. 


heights 


err , ' 


winds 


5.0 


10.0 


20.0 


0.5 




0.6842 


0.6824 


0.6925 






0.9328 


0.9307 


0.9332 






0.9473 


0.9469 


0.9479 


O 

H 




0.6843 


0.6824 


0.6924 






0.9329 


0.9306 


0.9331 






0.9473 


0.9468 


0.9478 


2.0 




0.6852 


0.6828 


0.6928 






0.9338 


0.9310 


0.9334 






0.9477 


0.9471 


0.9481 



Table 6. 3. 5-5: Variation of expected error 



Table 6. 3. 4-4: Variation of expected error 
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APPENDIX 2: 



SUBROUTINE TCOR 



This appendix contains the listing for the subroutine that 
evaluates the true covariance and cross-covariance function 
values between a given (grid) point and the observation points. 

— — — -f- — — * — — 4- — 

SUBROUTINE TCOR ( NCF , THD , PHD , NH , TUD , PUD , NU , TVD , PVD , NV , TD , PD , KH , KU , 

1 KV , VAR , CAPPA , CVM , NCVM ) 

C 

C THIS SUBROUTINE CALCULATES THE TRUE COVARIANCE BETWEEN THE 

C OBSERVATION VARIABLES AND INDICATED VARIABLES AT A GIVEN POINT 

C 

C THE ARGUMENTS ARE 

C 

C INPUT ARGUMENTS 

C NCF 

C THD , PHD ; NH 

C TUD,PUD;NU 

C TVD,PVD;NV 

C 
C 

C KH 

C KU 

C KV 

C VAR 

C CAPPA 

C TD,PD 

C NCVM 

C 

C OUTPUT ARGUMENT 

C CVM - ARRAY OF COVARIANCES BETWEEN OBS VARS AND GRID LOGS 

C AN NH+NU+NV BY 3 ARRAY 

C 

IMPLICIT REAL* 8 (A-H,0-Z) 

PARAMETER (NSZ=36) 

DIMENSION THD(NH) ,PHD(NH) ,TUD(NU) ,PUD(NU) ,TVD(NV) ,PVD(NV) , 

1 CVM(NCVM,3) 

DIMENSION TOH(NSZ) ,POH(NSZ) ,TOU(NSZ) ,POU(NSZ) ,TOV(NSZ) ,POV(NSZ) 
DIMENSION AP(5),AAP(5) 

COMMON /lO/KOUT 
DATA NTIME/0/ 

DATA AP , AAP/5D-6 , 3 . 0825D-6 , 3 . 0825D-6 , 3 . 0825D-6 , 2 . 188D-6 , 

1 ODO,ODO, .15D0, .2722D0, .2722D0/ 

C 

C THIS ROUTINE IS SET UP TO ACCEPT GENERAL (ISOTROPIC) COVARIANCE 

C FUNCTIONS. THIS REQUIRES THE DEFINITION OF A NUMBER OF ARITHMETIC 

C STATEMENT FUNCTIONS WHICH DEFINE THE DISTANCE FUNCTION IN TERMS OF 

C LATITUDE AND LONGITUDE (PHI AND THETA), ITS DERIVATIVES WRT PHI AN 

C THETA, AS WELL AS THE COVARIANCE FUNCTION AND IT'S DERIVATIVES. 

C 

C THE FUNCTIONS NEEDED ARE 

C DIST - DISTANCE FUNCTION 



COVARIANCE FUNCTION NUMBER 

ARRAY OF HEIGHT OBSERVATION LOCATIONS, NH OF THEM 

ARRAY OF U-WIND OBSERVATION LOCATIONS, NU OF THEM 

ARRAY OF V-WIND OBSERVATION LOCATIONS, NV OF THEM 

THE ABOVE ARE IN DEGREES, LONGITUDE AND LATITUDE, 
RESPECTIVELY 

NONZERO IF H COVARIANCE TO BE COMPUTED 

NONZERO IF U COVARIANCE TO BE COMPUTED 

NONZERO IF V COVARIANCE TO BE COMPUTED 

VARIANCE OF HEIGHT-HEIGHT ERRORS 
CORIOLIS CONSTANT 
GIVEN (GRID) LOCATION, DEGREES 
ROW DIMENSION OF CVM ARRAY 
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cov 

DCOV 

D2COV 

vAun 

VAUV 
i).‘:i)i> 
DSDT 
1)1)1' I 
DDI'Z 
1 ) 1 ) 1 ') 



ir.OTKOPrC COVAIMANCP, FUNCTrON AS A FUNCTION OF DISTANCE 

m;i<i VATivF OF COV win’ dtstanci': 

2ND DKRlVATrVE OF COV WRT DISTANCE - DCOV/DTSTANCE 

VARIANCE OF U-WIND 

VARIANCII OF V-WIND 

RARTIAI. DTST WRT I 'I IT I 

RARTFAI. DIST WRT TIIETAT 

PART OF 2ND PARTIAI. DIST WRT PIIII AND PIIIJ 

PART Ol’ 2ND PARTIAL DIST WRT TIIETAT AND TIIETAJ 

PAET of 2ND PARTIAL DIST WRT PUTT AND TIIETAJ 



DATA RAD/L. 372Dr,/ 

I ) I : ;t ( P , P I , T , 'P 1 ) - r ad* DSQRT ( ( P- P I ) * * 2 I ( ( T-T1 ) *COS ( ( P l pi ) /2 . do ) ) * * 2 ) 
COV(S) - VAR*(OMAA*( I .DO(A*S)*EXP(-A*S) (AA) 

DCOV(f:) ' -VAR*OMAA*A*A*S*EXP(-A*S) 

D2COV(r.) • VAR*OMAA*A**T*S*EXP(-A*S) 

VARU(P) - VAR*OMAA*(RAD*A)**2 
VARV(P) - VAR*OMAA*(RAD*A*COS(P) )**2 

I)SI)P(P,P1 ,T,TI ,S) <■' RAD**2/r.*(-(T-TI)**2*STN(PlPl)//l.D0l-(P-Pl) ) 

d;;dt(p,pi ,t/pi ,r.) - rad**2/s*(t-ti )*cos( (pipi)/2.do)**2 

DDPI ( P,PI ,T,TI ,S) « RAI)**2/S*(I .DO t ( T-T1 ) * * 2*COS ( P I PI ) /4 . DO ) 
DI)P2(P,PI ,T,TI ,S) - (RAD*COS( (PIPI )/2.D0) )**2/S 
l)l)P.)(P,PI ,T,TI ,S) - RAD**2/S*(T-T1 )*SIN(P)P3 )/2.D0 
I F(NTIME.NE.NCF) Til PIN 

flET PARAMETI'IRS FOR TUTS ITINCTTON NUMDER 



NCFR MAX(I,NCF) 
NCFR-MI N(NCFR, S) 

A - AP(NCFR) 

AA ■ AAP(NCFR) 

OMAA I . DO-AA 
WRITE (KOUT, I )A,AA 
PI -DATAN( 1 .1)0) *4. DO 
DTR - PI /I MO. DO 
NTIME - NCF 
END I F 



CONVERT LOCATTONf. TO RADIANS 

Dt) 100 1-1, Nil 

POII(I) - P1ID( I ) *I)TR 
TOIlfl ) TIID^l ) *I)TR 
100 CONTINUE 

DO I I 0 T-I ,NU 

POU( I ) ' PUI)(I)*I)TR 

TOU^I) • TUI)(T)*DTR 
MO CONTINUE 

DO 120 T«1,NV 

POV(I) « PVD(I)*DTR 
TOV(I) - TVDri)*DTR 
12.0 CONTINUE 

P PD*I)TR 
T • TD*DTR 



()0 



C START CALCULATIONS 

K = 0 

SP = SIN(P) 

CP = COS(P) 

IF(NH.NE.O) THEN 
DO 200 1=1, NH 
K = K + 1 

S = DIST(P,POH(I) ,T,TOH(I) ) 

DF = DCOV(S) 

IF(KH.GT.O) CVM(K,1) = COV(S) 

CVM(K,2) = O.DO 
CVM(K,3) = O.DO 
IF(KU.GT.O .AND. S.GT.O.DO) 

1 CVM(K,2) = DF*DSDP(P,POH(I) ,T,TOH(I) ,S) 

IF(KV.GT.O .AND. S.GT.O.DO) 

1 CVM(K,3) = DF*DSDT(P,POH(I) ,T,TOH(I) ,S) 

CVM(K,2) = -CVM(K,2)*CAPPA/SP 
CVM(K,3) = CVM(K,3)*CAPPA/CP/SP 
200 CONTINUE 
END IF 

IF(NU.NE.O) THEN 
DO 220 1=1, NU 
K = K + 1 

S = DIST(P,POU(I) ,T,TOU(I) ) 

DF = DCOV(S) 

DD2 = D2C0V(S) 

SPI = SIN(POU(I) ) 

CPI = COS(POU(I)) 

CVM(K,1) = O.DO 
IF(KH.GT.O .AND. S.GT.O.DO) 

1 CVM(K,1) = DF*DSDP(POU(I) ,P,TOU(I) ,T,S) 

CVM(K,2) = VARU(P) 

IF(KU.GT.O .AND. S.GT.O.DO) 

1 CVM(K,2) = -DF*DDP1(P0U(I) ,P,TOU(I) ,T,S) + 

2 DD2*DSDP(P,P0U(I) ,T,TOU(I) , S ) *DSDP ( POU( I ) ,P,TOU(I) ,T,S) 
CVM(K,3) = O.DO 

IF(KV.GT.O .AND. S.GT.O.DO) 

1 CVM(K,3) = -DF*DDP3(P,P0U(I) ,T,TOU(I) ,S) + 

2 DD2*DSDP(P0U(I) ,P,TOU(I) , T , S ) *DSDT( P , POU( I ) ,T,TOU(I) ,S) 
CVM(K,1) = -CVM(K,1)*CAPPA/SPI 

CVM(K,2) = CVM(K,2)*CAPPA**2/SP/SPI 
CVM(K,3) = -CVM(K,3)*CAPPA**2/SP/CP/SPI 
220 CONTINUE 
END IF 

IF(NV.GT.O) THEN 
DO 240 1=1, NV 
K = K + 1 

S = DIST(P,POV(I) ,T,TOV(I) ) 

DF = DCOV(S) 

DD2 = D2C0V(S) 

SPI = SIN(POV(I) ) 

CPI = COS(POV(I) ) 

CVM(K,1) = O.DO 
IF(KH.GT.O .AND. S.GT.O.DO) 
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240 

1 



CVM(K,1) = DF*DSDT(POV(I) ,P,TOV(I) ,T,S) 

CVM(K,2) = O.DO 
IF(KU.GT.O .AND. S.GT.O.DO) 

CVM(K,2) = -DF*DDP3(POV(I) ,P,TOV(I) ,T,S) + 
DD2*DSDP(P,POV(I) ,T,TOV(I) , S ) *DSDT ( POV( I ) ,P,TOV(I) ,T,S) 
CVM(K,3) = VARV(P) 

IF(KV.GT.O .AND. S.GT.O.DO) 

CVM(K,3) = -DF*DDP2(P,P0V(I) ,T,TOV(I) ,S) + 
DD2*DSDT(POV(I) ,P,TOV(I) ,T , S ) *DSDT( P , POV( I ) , T ,TOV( I ) , S ) 



CVM(K,1) 

CVM(K,2) 

CVM(K,3) 

CONTINUE 

ENDIF 

RETURN 

FORMAT (' TCOR 
END 



CVM( K , 1 ) *CAPPA/CPI/SPI 
-CVM( K, 2 ) *CAPPA**2/CPI/SPI/SP 
CVM( K , 3 ) *CAPPA**2/CPI/CP/SP/SPI 



- 01/12/90: 2ND AR, A,AA = ',1P,2E11.3) 
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